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I. GENERAL INTRODUCTION 


This document constitutes the First Quarterly Report submitted 
by ADCOM, Inc. to the NASA Goddard Space Flight Center under Contract 
No. NAS 5-9011. The work covered here was performed at ADCOM, Inc. 
between 1 July 1964 and 30 September 1964. 

The primary objective of this program is to demonstrate the feasi- 
bility of Reducing the noise threshold for demodulating FM signals at least 1 
3db below the threshold experienced with a second-order phase-locked loop 1 
demodulator whose design is optimized for best threshold performance in 
demodulating a given signal. 

The technical plan of the program calls for a critical assessment 
of the current state of knowledge of the most promising threshold reduction 
approaches, taking into full account the available technical literature. The 
object of this assessment is to identify the various mechanisms that cause 
the threshold effect and limit the ultimate threshold reduction accomplish- 
ment with each of the various promising demodulation techniques. Inasmuch 
as the literature on this subject abounds in widely quoted analyses and re- 
sults of a controversial nature, critical reviews of some of these publica- 
tions are in order. 
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In the light of the results of the above critical examination of the 
various approaches, a few of the most promising approaches are to be 
singled out for implementation and thorough testing. 

Section II of this report presents a synopsis of results reported 
in the main body of the report. 

Section III presents a critical treatment of certain fundamental 
aspects of conventional demodulators, phase-locked loop, frequency 
compressive feedback, and band-dividing demodulators that are essen- 
tial in the analysis of performance, evaluation of improvement poten- 
tialities, and comparison of ultimate performance limitations of each of 
these techniques. 

Section IY presents a critical review and evaluation of the validity 
and significance of certain published analyses, and an extension of the 
results of one of the reviewed papers. 

In Section V, analyses are presented of higher-order phase-locked 
loops to determine their ultimate optimized performance capabilities and 
requirements. 

In Section VI, some aspects of experimental test and evaluation are 
discussed and an automated procedure is described. 

The work plans for the next interval are outlined briefly in Section VII. 
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II. SYNOPSIS OF REPORT 

The present chapter summarizes the main investigations reported 
in the sections that follow. 

Chapter III characterizes the basic models of the PEL and FCF 
demodulators. The PLL is first characterized as a nonlinear, randomly 
time-variant system and the necessary assumptions involved in the con- 
ventional linear time -invariant model are indicated. An analogous analysis 
is then effected for the FCF demodulator. The absolute lower limit on 
closed-loop noise bandwidth is established for both systems and threshold 
reduction capabilities are considered on this basis. Band-dividing demod- 
ulation is finally considered and the limitations of different signal proces- 
sing techniques are established. The case of CW signals is shown to yield 
a bandwidth comparable to that of a second-order PLL demodulator for each 
of the filters in the band-dividing bank so that the SNR performance of the 
Akima system will not show the desirable improvements. The case of 
sampling and quantizing is also eliminated because of the large number of 
filters required leaving sample-and-hold processing as the other possi- 
bility to be considered. 

Chapter IV evaluates certain published analyses of interest. A 
detailed discussion of some of the limitations of Develet's papers is first 
presented. The work of Sanneman and Rowbotham is then evaluated and 
some basic results are clarified and extended. 
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Chapter V evaluates the performance of a second-order PLL 
demodulator and analyzes the possible threshold improvements attainable 
with higher order PLL's. The second-order loop performance is first 
optimized for a sinusoidal modulation and the above -threshold operational 
requirement is established. An analogous approach is used for the third 
order loop and the improvement capabilities are identified and found to 
be limited from feasibility considerations. The optimum realizable PLL 
demodulator is then determined for' the cases of rectangular and an RC 
PM spectra. The threshold characteristics of these receivers are 
established and possible approximations to the optimum transfer function 
are finally presented and evaluated. 

Chapter VI is concerned with the experimental program. The 
previous A DCOM work of SNR measurements is summarized and a 
threshold test system for the present contract is suggested. The nec- 
essary modifications that will yield an automatic, rather than manual, 
threshold curve plotter are included. 
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HI. BASIC MODELS OF EXPONENT DEMODULATORS 

In the analysis and comparison of threshold reduction techniques, 
it is of the essence that the basic models and assumptions used correspond 
to the physical conditions encountered in the real world. In. the last few 
years, a number of attempts have been made at determining ;and co m par in g 
the ultimate performance of various exponent demodulation techniques. Un- 
fortunately, many of these publications have introduced and perpetuated a 
host of now widespread misconceptions about thebasic theoretical models and 
the fundamental considerations that define the performance characteristics, 
limitations, and comparative merits of these techniques. 

Inasmuch as the major objective of this program is to demonstrate 
the feasibility of reducing the noise threshold of exponent demodulation at 
least three db below the threshold experienced with a second- order phase- 
locked loop, it is important that we re-examine the basic models and con- 
siderations that are essential for predicting the ultimate threshold reduction 
performance of modulators. This re-examination is essential for dispelling 
certain published analyses from the picture lest their results be considered 
to define true limits inhibiting further consideration of otherwise promising' 
approaches. On the positive side, the re-examination brings into sharp focus 
the basic criteria and causes of the noise thresholds of various demodulators. 
In this way improved insight into the avenues for improving performance 
should result. 
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3. 1 The Conventional FM Demodulator 


The conventional limiter-discriminator demodulation technique 
provides a convenient standard of comparison for other demodulators. The 
performance characteristics of this type of demodulator have heen the subject 
of many analyses and discussions, and hence will not be taken up here. How- 
ever, one point -is considered to be of sufficient importance here to deserve 
presentation. 

The point of interest concerns the common assumption about the de- 
pendence of the threshold of conventional FM demodulation upon the deviation 
ratio of the anticipated FM signal. Thus it is usually correctly stated that 

the threshold S/N of the FM demodulator varies directly with the noise band- 
width of the i-f filter that shapes the input noise spectrum. However, it is 
unjustly stated that the (noise) bandwidth of the i-f filter varies directly with 
the frequency deviation of the expected signal, and hence the noise threshold 
of the conventional FM demodulator varies ultimately directly with the de- 
viation ratio of the expected FM signal. 

Actually, it can be shown in a number of ways that the bandwidth of 

a given type of filter varies with the square-root of the frequency deviation 

of the expected signal, if the tolerance on distortion and the type of filter are 

1 

held fixed. In fact, direct analysis in terms of series expansion of the re- 
sponse of a linear filter to an FM excitation shows that for a given type of 
filter and a specified tolerance on acceptable distortion the filter bandwidth 
(BW) is given by 
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(BW> = K */ | <i. (t) j (3.1) 

' i ‘max 

where uk (t) = dw^(t)/dt, and wdt) is the instantaneous frequency of the input 
FM excitation. The constant K is determined by the filter type and the tol- 
erance on distortion. Thus, if 


w..(t) = A£2g(t), jg(t)j 


max 


= 1 


(3.2) 


then 


CO . (t) I = Aft|g(t)| 
l 'max 1 & ‘max 


and 

(BW) - (3.3) 

The proportionality (3. 3) can also be shown as follows. We reason 
that as the FM signal frequency sweeps across the filter passband, the 
filter response will depend upon the sweep rate. In particular, if the sweep 
rate is such that the signal lingers within the nominal passband of the filter 
for a length of time equal to k time constants of the filter, then if k > 2 or 
3, the filter response should build up to anear-steady-state pattern, the desired 
closeness of the pattern to the steady-state form dictating how high k should 
be. But 

(filter time constant) = k^/CBW), k = const. 

and 


(maximum sweep rate) = k 2 AO , k^ = const. 

Consequently, for a specified quality of quasi- static response, we must have 

(BW)/k AU = signal dwell time within the nominal passband of the 
filter 

= kk^BW) 
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whence 

(BW) « ^AO , Q.E.D. 

It is a simple numerical exercise to show that use of the commonly 
employed formula 

(BW) = 2(An+co m ) (3.4) 

where co^ = highest modulation frequency, results in a level of modulation 
distortion that varies with A £2, for a fixed type of filter. 

It has been argued that using the narrower bandwidth indicated by (3.3) 
would cause the FM signal to partially sweep the falling skirts of the filter 
response, thus causing its amplitude to dip relative to the noise over parts 
of the modulation cycle. But the important fact is that the total noise passed 
by the narrower filter is proportionately lower than the total noise passed by 
a filter designed in accordance with Eq. (3.4), so that a net reduction in the 
noise threshold would be perceived in a conventional FM demodulator . 

In fact, the point may be stretched farther to suggest that one may 
consider designing the i-f filter so that its bandwidth is reduced at very low 
signal levels so that distortion in the filter is traded for a lower demodu- 
lator noise threshold. In many applications signal distortion is more pal- 
atable than the considerably sharper degradation encountered below the 
demodulator noise threshold. 

The above observations are of fundamental importance in the evaluation 
and comparison of the threshold performance of various FM demodulators rela- 
tive to the performance of a conventional FM demodulator. 
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3. 2 The PLL Demodulator 


In the study of the ultimate capabilities of PLL demodulators, two 
questions are of paramount interest in this report. The first is the absolute 
lower limit on the value of closed-loop bandwidth, the second is the modeling 
of PLL demodulators for analytic purposes. Both of these questions will be 
reviewed here for support of later arguments. 

3.2.1 Basic Models ^ 

The basic functional block diagram of a phase -locked loop 
(PLL) demodulator is shown in Fig. 1. The signal plus noise at the output 
of the i-f amplifier may be expressed as 

e. f (t) = [E g + x^ ^(t)] cos [co^t + $(t)] 

" X q, jf^ sin + $ W] (3. 5) 

in which 

E g = constant signal amplitude 


if) (t) = instantaneous phase fluctuations of 
desired signal, 

« (message time function) in <{>M, 

cc (integral of message time function) 
in FM, 


and the noise is represented as the sum of two components -that are instantan- 
eously cophasal and in quadrature with the modulated signal and have instan- 
taneous amplitudes x (t) and x . f (t), respectively. The VFO oscillation 

0 j Xjl q j ll 


may be expressed as 


e (t) = E cos [<o t + <J) (t)] 

osc osc osc T osc 


(3.6) 
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e in (0 


Fig, 1 




Basic functional block diagram of 
phase -locked loop demodulator. 
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The oscillator instantaneous phase fluctuations, <b (t), in general combine 

osc 

fluctuations due to desired signal modulation, attendant noise and nonlinear 

interactions of signal and noise. Under the assumption of stable feedback 

operation, and in the absence of spurious' by-products in the output of the 

multiplier, we have, with co. = w 

if osc 

^ iosc (t) = {t 1 + i Lx c,tf (t) ] Si " [* (t) - hsc (t) ] 
r s 

+ V if (t) COS [* (t) ' +osc (t) ]} H h lp (t) <S - 11 
s 

where $3 denotes convolution and (t) is the impulse response of the low- 
pass filter. The existence of a feedback steady-state condition in which 
Eq, (3. 7) gives an adequate description of loop performance requires that 

| A (j) ] = j lb (t) - <j) (t) j < 7r/2 (3.8) 

Equation (3. 7), subject to condition (3. 8), suggests the model shown 
in Fig. 2 for PLL operation. If j A<|> ] « 1 almost all of the time, then the 
sine term may be approximated by its argument and the cosine term by unity. 

In many practical situations, it is sufficient to require [ A <j) | <[ ~ 6 almost 

all of the time in order to permit the approximation sin A t|> « A<j>, and | A<j>| 

< y 0. 2 almost all of the time, in order to permit the approximation cos A <|> 

& 1, In such situations, the distortion ignored by these approximations is 
considered to fall below tolerable limits. (This, however, is not true in 
general, especially in demodulating multichannel FDM/FM signals. ) 

It is important to note that the condition for approximating cos A <j> by 
unity is more stringent than that of approximating sin A<j> by A<j) . Consequently, 

11 
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Fig. 2 Model of PLL demodulator in the absence of pre -limiting of 
sum of signal and noise, and under the assumptions of stabli 
loop operation, and lifa(t) - J> (t) I < ir/2. 







AD, 


CO 


if cos A<j> is approximated by unity, then it follows necessarily that sin A<j> 
Aij) and A cannot be allowed to exceed */ 0. 2 too frequently. 

Thus, approximation of the model of Fig. 2 by the linear model of 
Fig. 3 requires that the following conditions be satisfied: 

(a) The factor at point^T) in Fig. 2 is approximately 

equal to unity. This means that the cophasal com- 
ponent of the i-f noise must be negligible, which 


is practically the case if 
(S/N) < 5 (or 7 db> 

sin [0(t) - (t) ]= sin A <j> ^ A<1> 


(3. 9) 


cos A <j> ^ 1 

These approximations hold if 

| A<j> ] « 1 almost all of the time (specifically, (3. 10) 

j A <f> | _<: V 0. 2) The condition | A.<}) | « 1 reduces 


the time function at point in Fig. 2 to i/>(t)-({) (t) 

and that at point (?) to x„ .„(t)/E . 

q, if s 

When the linear model holds, the output of the low -pass filter is given 
by the inverse transform of 


/iH (s) 

V- - s *osc is)la - 3 + M aH ip (,) L s * (s> + E ; s Vf (s) ] (3 - n) 
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Fig. 3 Baseband linear model of phase -locked 
loop demodulator operating upon sum of 
FM signal plus .noise. It holds only sub- 
ject to the conditions 

i^ (t) - <j> (t)j « 1, <«> 

OS c , s 

and (input S/N in the i-f bandwidth) > 7 db. 
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Thus the use of a low -pass filter that does not cut off abruptly enables fre- 
quencies outside the minimum necessary band to be incompletely attenuated, 
thereby widening not only the equivalent noise bandwidth of the low -pass filter 
itself, but also the equivalent closed-loop bandwidth for the transmission of 
instantaneous frequency or phase perturbations from the input signal termi- 
nal to the VFO. Low -distortion PLL operation requires that the desired base- 
band waveform j/)(t) be transmitted to the instantaneous frequency of the VFO 
with negligible distortion. But this does not necessarily mean that the loop 
low -pass filter need pass the desired modulation, $(t), with negligible dis- 
tortion because, within the loop, this filter operates on 

jus ^ (s) 


s + JUaH (s) 


which may in general be a linearly predistorted form of s#(s), but the feed- 
back factor [i a must be such that (s) will provide proper linear compensa- 
tion so that the resultant system function 

tlaH.Js) 


JJL 


H eq (s) s + juaH ip (s) 


(3. 12) 


will yield 0(t) with negligible distortion. Under these conditions, Eq. (3. 11) 
can be simplified to the form 


-T, s 


s$ (s) ~s*(s)e d,eq +^-sX .is) H (s) 


osc 


-L-l- \ 

g eq 


in which r , ^ is the time delay in distortionless transmission through H (s). 
n, eq eq 

Inverse transformation yields 
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• (t) 
r osc 




4~ X .it) g? h (t) 

E q, if eq 


(3. 13) 


This shows that under conditions of high input S/N ratio, the phase -locked 
loop simulates the bandpass analog of its equivalent closed-loop instantan- 
eous phase or frequency low -pass transfer function, with its center fre- 
quency automatically and continuously adjusted to follow closely the instan- 
taneous frequency of the input signal, provided that the closed-loop equivalent 
low -pass filter can pass the signal modulation function with negligible dis- 
tortion and conditions (3.9) and (3. 10) are satisfied. 

It is important to point out that as long as condition (3. 10) is satisfied 
to within acceptable tolerances, the condition that H^(s) pass the desired 
modulation waveform if) (t) with negligible distortion is not necessary for 


proper FM demodulation. In fact "optimization 11 of loop performance with 
} 

•respect to both noise and distortion above the threshold of linear variation of 
output S/N ratio with i-f S/N ratio may result in non -negligible linear dis - 
tortion of i/)(t) by H^{s). Such distortion can, however, be reduced to within 
tolerable bounds by linear compensation, or "equalization, 11 of the offending 


non -uniformities in the characteristics of H (s) by means of an appropriate 
i eq 

low-pass filter after the loop. Such a filter is also generally desirable from 
the viewpoint of baseband noise filtering because H^/s) will usually have a 
|arger noise bandwidth than is necessary in proper overall low -pass filtering 
of the baseband waveform. 


I 
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If a perfect bandpass limiter operates on the sum of signal 
and noise before it is applied to the APC loop, Eq. (3. 7) is replaced by 


— <b (t) = h (t) sin U<t) + 9 .it) - <1> (t) T (3. 14) 

jua r osc ip L r n, if y osc J 

x .it) 


where 


0 .M = tan' 1 - * - f — v; 
n ' 11 E s + X c, if (t) 


(3. 15) 


This suggests the model shown in Fig. 4(a). Approximation of this model by 
a linear model requires that 

|0(t) + 0 n ^.(t) - I «1 almost all of the time. (3.16) 


Moreover, if only the quadrature component of the i-f noise is to be included 


in the input, then we must also have 

9 .it) ^T“X .it) 

n, if E q, if 


(3. 17) 


which holds provided that (S/N) exceeds or equals a "practical" threshold of 10 db. 
Alternately, Eq. (3. 14) maybe written as 

— i (t) = h„ (t) sin [0(t) - <{> (t)l cos 0 .it) 

Hot y osc ip I ^ r osc n,it 


+ cos [ ip(t) - <{) (t)] sinS .it)} (3. 18) 

r T osc n,if J 

which suggests the model in Fig. 4(b). Again, linearization requires that 
condition j A<j) | « 1 and (S/N)^ > 10 db. 
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Fig. 4 Models of PLL demodulator when the input is the amplitude - 
limited resultant of signal and noise , under the assumptions 
of stable loop operation and | j/>(t) - £ (t) \ £ W 2. 

OSC 
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3 - 2. 2 Absolute Lower Limit on Closed-Loop Noise Bandwidth 

Let us now consider the situation in which the loop lowpass 
filter is assumed to have uniform transmission over B (« ua) rad/sec-- 
the frequency range occupied by the baseband modulation waveform g(t) -- 


with a constant delay sec and with very large attenuation outside the 
stated bandwidth. It is well-known that such an attenuation characteristic 
cannot be realized physically without a correspondingly large delay that is 
approximately constant only over a fraction of the passband near the center 
frequency but rises rapidly near the edges of the band. Such a large delay 
can be shown to be totally unacceptable for realizing a stable feedback sys- 
tem in which the PLL can lock to an FM signal, and yield a negligibly dis- 
torted replica of the modulation function in its output. But the idealizations 
we .make concerning constant low delays associated with .very sharply selec- 
tive attenuation characteristics are merely intended to be a crude imitation 
of physical systems that helps in determining the absolute lower bound on 
the closed-loop noise bandwidth. 

Thus, with the input excitation described by Eq. (3. 5) and 
the oscillation expressed as in Eq. (S.'S), then, in the absence of spurious 
by-products in the output of the multiplier (or "phase detector"), the low- 
frequency components applied to the lowpass filter are 

r x c if ^ q i 

+ — E J Sin ^ + e" 4 x „ Jfft) GOS a< !> (3.19) 
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A feedback steady-state condition exists, in which (3, 19) properly describes 
the low-frequency components applied to the lowpass filter, provided that 


| A<J>| < W 2. If we assume that Acj> < J 0. 2 , then (3, 9) can be simplified to 
the form 


4 1+ ^- X c,if (t) ] ^ x q,if (t> 

Cl c-1 ± 


The output of the lowpass filter will therefore be 


e ip (t) „ M [*<t - T ip ) - ^ sc (t - T^)] 


(3.20) 


+ # K, if (t> [ - ■t > oso (t) ] }® V t,+ #■ X q,if (t) B V‘> 


ip' 
(3.21) 


If the baseband bandwidth is a small fraction of the i-f bandwidth, the second 
term in Eq. (3. 21) can be approximated by 


- *osc (t) ] [ x c,if (t! B h ip (t> ] 


(3. 22) 


When expressed in this form, it becomes clear that the second term in Eq. 
(3. 21) is practically negligible if the ratio of mean squared signal to mean 
squared noise contained within twice the idealized lowpass filter bandwidth, 
centered about the instantaneous frequency of the input signal, exceeds or 
equals '5( db). Expression (3. 20) shows that, in general,, the second term 
in Eq. (3. 21) can for practical purposes be neglected as long as the i-f S/N 
ratio exceeds or equals 7 db. 
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If the second term in Eq. (3. 21) can be neglected, we have 

%(t) “ “ O - V - <>osc (t - Vl + W~ x q ,if (t > B %<« ,S - 231 

Thus, under high S/N conditions we have 


= a % (t> 

83 - V - ^ osc (t - V + IP \ If (t> B V (t> < 3 - 24> 

s 


whence 

s^(s) « *( B ) e' s/ ^ a + i s x (s) e" s/Ma . H (s) 

E s £p 


(t) M A Og(t - 1 1 (icy.) + X (t - 1 /m a ) ® h (t + T ) 
ot)C q,ll i!p x. p 

s 


(3. 25) 


This shows that the delay required for yielding FM demodulation by the 
phase-locked loop with a sharp cutoff lowpass filter is given byT/jj a, the 
time constant of the loop excluding the normalized, idealized lowpass filter. 


In order to guarantee that this delay is small but nonzero so that proper FM 
demodulation can result, the loop gain fj,a should be large (but not infinite). 

The result expressed in Eq. (3. 25) shows that the phase- 
locked loop with a rectangular lowpass filter within the loop has a noise 
bandwidth equal to the bandwidth of the lowpass filter. But the minimum 
permissible value of this bandwidth is twice B^(in order to include positive 
and negative frequencies) where is the (positive -frequency) bandwidth 
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of the desired message. Consequently* we may conclude that the absolute 
lower bound on closed-loop noise bandwidth is 233^* and* hence* the abso- 
lute upper bound on the threshold reduction factor achievable with a phase - 
locked loop is 


2* (B /2B ff ) = B.„/B 
if Ip if ip 


(3. 26) 


The factor of 2 in (3. 26) takes into account the fact the threshold of linear 
variation of baseband S/N ratio with i-f S/N ratio for a PLL> demodulator is 
reached when the S/N ratio in the oscillation noise bandwidth equals 5(^7 db). 
3. 3 The FCF Demodulator 

The functional diagram of an FCF demodulator is shown in Fig. 5. 

We assume that the closed-loop system can lock properly to the desired signal 

in the absence of noise, and let e.„(t) and e (t) be described by (3. 5) and 

it osc 

(3. 6)* respectively. The FCF loop is assumed to be driven with no prior 
amplitude limiting of signal plus noise. Under these conditions, we further 
assume that the mixer delivers in the nominal passband of the loop bandpass 
filter only the beat components 

[ E S + X C, if (t)] 003 [ “bp* + * (t> - *030^ 

£ q, if(*> Sin [w bp t + < ,(t) '- , l’o S c (tn 


- X 


in which co = co. - w = nominal center frequency of the bandpass filter 
bp if osc 

within the loop. If an amplitude-insensitive linear FM demodulator and a 
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Fig. 5 Basic functional block diagram of frequency- compressive feedback demodulator. 
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linear feedback from the discriminator output to the instantaneous frequency 
of the VFO are assumed, we can show that 


<j> (t) 


osc 



A q (t) EJ h^t) 

A (t) E5 h T _<t) 
c DP 



(3.27) 


where 

A c ( t) - [l + 5 - x 0 jif <t>>os [0 (t) - 4> <t)] 

s 

- X q, if (t) Sin Asc™] 

s 

A q (t) ' [ 1 + IT X o,lt (t) ] Sln [* (f) - +osC (t> ] 

s 

+ x q, if (t) cos [ 5 |)(t) - ^ osc (t) ] 

and h T (t) is the impulse response of the low -pass analog of the bandpass 

J-iJr ^ 

filter. 

3. 3. 1 Basic Models 2 

Expression (3. 27) for <1* can be used to derive lowpass 

equivalent models of the FCP operation upon the input signal plus noise for 
arbitrary input SNR’s. We illustrate this with a number of important exam- 
ples. 

Consider first the situation in which the noise is so weak 
that it can be ignored entirely (i. e. , the noiseless case). Here (3. 27) is 
simplified to 
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<i> (t) 

osc 


A (t)Bsin[i|)(t)-|) (t)] 

k /? h. (t) 33 tan r- - ^ i - , 

d K fb h ._(t) ®cos[)/)(t) - <b (t)l 

Lif • osc 


This suggests a general nonlinear baseband model of the FCF operation upon 
signal modulation. It is clear, however, that linear-approximation models 
of the FCF operation upon the signal modulation can be justified under either 
of the two conditions: 


(a) 


U(t) - <J> os (t) I « 1 (or < 0. 2 for many practical 


purposes) which allows the approximation 


tan 


-1 ^LP^ 33 sin [^>(t) - 4> ogc (t)] 
h (t) EJ cos [ ih(t) - <{) (tj] 

L-iir OSC 


\ P (t) s [#(tl - * 0BC (t)] 


( 3 . 29 ) 


(b) 


or 

the response of the loop bandpass filter to the fre- 
quency-compressed signal is quasi -stationary with 
negligible nonlinear distortion of the modulation, 
which allows the approximation 


tan 


-1 h LF^ S Sin ” ^osc^ 

h Lp (t) 3S cos [fC t) - ?~(t)] 


osc 


tan 


-1 


sin[0(t - r d ) - 4> (t - r.) 1 


osc 


008 [* (t - r d ) - 1» oso (t - T d )] 


* (t ‘ T 1 - +oso (t ‘ 


( 3 . 30 ) 


Condition (a) leads to the model of Fig. 6(a); condition (b) 
leads to 6(b). The two models differ only in the equivalent lowpass analog 
of the loop bandpass filter 
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(a) Subject to condition that the FCF reduces the phase deviation 
at the input of the loop bandpass filter to a value much smaller 
than unity (0. 2 rad or less). 



(b) For operation on signal modulation. Valid only when the loop 
bandpass filter responds to the frequency- compressed signal 
in a quasi- stationary manner, with negligible distortion, and 
delivers replica of compressed modulation function delayed 
T, sec. 


Fig. 6 Linear models of FCF demodulator. 

26 

ADVANCED COMMUNICATIONS • RESEARCH AND DEVELOPMENT 









is given by 



With reference to Fig. 6, the closed-loop system function 


k hW s,H ip< s) 


FT (o’! — 
aq' - 1 + 


(3. 31) 


for the model of Fig. 6(a), and 


H (s) = 
eq 


— 7* S 

Vtb e di w s) 

1 + v® e TdS v s) 


(3. 32) 


for the model of Fig. 6(b). 

The frequency modulation of the signal at the input of the loop bandpass 
filter is given by the inverse transform of 


H (s) s^(s) 
eq, c 


(3. 33) 


in ■which 


H (s) » 1 - H (s) 
eq, c eq 


(3. 34) 


is the complement of the closed-loop system function H e ^(s). In general, 
ma y cause linear distortion of the desired modulation, but such 
distortion must not exceed certain low bounds if proper operation is to be 
maintained in the presence of noise down to and below the threshold of a 
conventional FM demodulator. 

When the linear models cannot be justified, the nonlinear 

effects of the FCF operation introduced in the response of the loop bandpass 

filter can be computed with the aid of Eq, (3. 27). 
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Linear and nonlinear approximation models of the FCF 
operation can also be developed when the noise cannot be ignored. Of par- 
ticular interest are the models that apply under conditions of high input S/' 
ratio. We distinguish the two situations: 


| ti>(t) + 9 .„(t) - <j> (t)| < 1 

n, if T osc 1 


almost all of the time and 




a considerable fraction of the time. 


In the first of these situations. 


TT^ *08c (t) - [ * » + 9 n,if (t) - *osc (, « Bh BP (t, ® i 1 p (t ' 

which establishes the applicability of the linear model of Fig. 6(a) with 
if) (t) + © n ^(t) as representative of input signal plus noise. This linear 
model yields a practically satisfactory approximation to the actual S/N 
mean-square ratio performance provided that 


6 n, if (t) ** E X q, if (t) 
s 


to within a first-order mean-square contribution from (1/EJ .„(t) that 

s c , it 

is at least 10 db below that of (l/E )x ,„(t) — i. e. , as long as the S/N 

s q,il 

ratio in the pre-FCF bandwidth exceeds or equals about 12 db. 
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For the situation in which (3. 36) applies, we recall that under 
conditions of high input S/N ratio, the effects of signal modulation and of noise 
upon the output can be computed separately. Thus, if we first ignore the noise, 
then with 

U<t) - <J> (t) I > 0. 2 

x osc 

the model of Fig. 6(b) applies subject to the condition that the loop bandpass 
filter response reproduce the compressed modulation with negligible distor- 
tion. Next, ignoring the signal modulation, we observe that with 


I V if tt) 1 " !'¥' X q,if (t) l 
s 

well below unity most of the time (i. e. , with the pre-FCF S/N ratio 

exceeding about 12 db), the linear model of Fig. 6(a) will approximate 

closely the effect of FCF upon the FM noise. Thus, we have one model for 

signal modulation(Fig. 6(b)) and another for noise (Fig. 6(a)). 

The preceding results show that the FCF demodulator 

operating upon the sum of signal plus noise simulates a linear tracking 

filter given by the bandpass analog of the equivalent closed-loop lowpass 

transfer function H (s) defined in Eq. (3. 31), with its center frequency 

eq 

continuously and automatically adjusted to tracklhe instantaneous frequency 
of the input signal, provided that the desired signal modulation is repro- 
duced with negligible distortion by the closed loop system, and (S/N)^. > 

12 db. 
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3. 3. 2 Absolute Lower Limit on Closed-Loop Noise Bandwidth. 

Let us now consider the situation in which the loop bandpass 

and lowpass filters in Fig. 5 are assumed to have uniform transmission 

over B_._ and B„ rad/sec, respectively, with constant delays of r , and 
BF Jcp Q 

* 1 - sec, and- with very large attenuation outside the stated bandwidths. Assume 
that relatively weak noise is superimposed on an input unmodulated carrier. 

Now, since 

a) the noise in the passband of the loop bandpass 
filter is much weaker than the signal; 

b) the open -loop lowpass filter action, by assump- 
tion, restricts the transmission to a bandwidth 
that is much narrower than the bandwidth of 
the receiver stages that drive the FCF demodu- 
lator; and 

c) the bandwidth of the loop bandpass filter is suf- 
ficient to enable this filter to respond quasi- 
stationarily to the weak frequency fluctuations 

by the noise (i. e. , Hpp(joj) is essentially uniform 
— with unit amplitude and phase slope say — 
over the range within which the noise-deviated 
instantaneous frequency of the local oscillation 
is bounded almost all of the time); 

we have at the output of the loop bandpass filter 

- [ E s + X C. rf (t > ® h LP (t> ] ' COS [“if* - hsc (t ‘ T d>] 

' [ X q.rf (t) Sh LP (t) ] ' Sin [“if* " 'W-V] 

(3. 38) 

In response to epp(t) , the amplitude -insensitive FM demodulator delivers 
to the loop lowpass filter an excitation proportional to 


=T x (t) & h (t) -A (t-r ) 
E q, rf LP osc d 

S n a 


(3.39) 
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If the loop lowpass filter has an impulse response denoted by h (t) , then 
k d 

e (t) = ■=— x It) i3 h (t) h. (t) - kj<j> (t - r - r ) 
ip E q, rf' LP ip d Y osc d ipf 

s 


TT“ n (t) - k <{> (t - T , ~ T ) 
E^ ip d T osc d ip 


s 


"3. 40) 


wherein we have defined 


% w -\rt ws V* ,B V) 


(3. 41) 


We have also chosen to ignore the filtering, except for delay r , of 
c|> (t) by h (t) in (3. 40) on the basis that <j> (t) is already a lowpass 

OSC i!p 


' OSC 


function that has received sufficient filter action by h„ (t) to allow it to 

J ip 

pass again through this filter with no modification. 

If we set (j> (t) = B„ e„ (t) in (3.40), we obtain 

osc fb ip 


V s) = 


k j/ E 

d ; s 


1+ Vtb e “ " p 


-<T J + T,_)S * N ip (S) 


which, for negligible open-loop delay and/or for (the tacitly assumed) 
large feedback, leads to 


V t) = TTF^ V t+T,) 


(3. 42) 


T> * 


k d^fb 

(Td + V 
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Finally., the noise transmitted to the oscillator phase is 


k /E 

<j> (t) = - - — | — n (t+r) 

osc 1 + k d ip 


(3. 43) 


From Eqs. (3. 41) and (3. 43) it is clear that the closed-loop noise band- 
width of an FCF demodulator with a rectangular lowpass filter within the 
loop is given by the bandwidth (including positive and negative frequencies) 
of the lowpass filter (i. e. , 2 33 ). This is then the absolute lower limit 

*-xr 

on FCF closed-loop noise bandwidth. 

The determination of the noise threshold requires an 

examination of which of two conditions will break down first: 

(i) Operation of the loop discriminator above 
its threshold; or 

(ii) (t) <1/20 (which guarantees that the 

osc., n 

contribution from the interaction of oscillator 

phase noise, (t (t) , and input noise is 
osc, n 

negligible). 

If the threshold defined by condition (i) is denoted (S/N) 1 
and that defined by (ii) is (S/N) 2 , then 

2 33 

(S/NF/tS/N) = < 1 (3.44) 

1 2 l^p 


which shows that condition (i) breaks down first and, hence, defines the 
closed-loop threshold. Consequently, the absolute upper bound on the 
threshold reduction factor achievable with FCF is 

B lf ' b bp (3 - 45) 
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In comparison with the upper bound for PLL, the upper bound given for 
FCF by (3. 45) is lower by a factor of 

B Bp /% >2 < 3 ' 46 > 

The limit of threshold reduction is achieved with FCF 
using idealized rectangular filters when the frequency compression reduces 
the necessary bandpass filter bandwidth to a value close to twice the base- 
band bandwidth. The upper bound on the permissible value of and 

the gradual cutoff encountered (and required) in practical filters restrict 
the degree to which this limit can be approached. 

3. 4 Feedforward Band-Dividing Demodulator 

One promising band-dividing scheme is proposed in Fig. 7. The 
instantaneous frequency of the input signal is estimated by a bank of filters 
or its time-domain equivalent. An FM signal at the estimated instantaneous 
frequency is then generated. The first mixer takes the frequency difference 
between the input signal and the estimated signal. If the estimate is reason- 
ably accurate (this point will be expanded later) then the frequency difference 
will have much smaller excursions than the input frequency. The difference 
signal may therefore be passed through a relatively narrowband filter, thereby 
eliminating much of the input noise power. The filtered difference signal is 
mixed with the estimate a second time, restoring the original excursions in 
input instantaneous frequency. The result of these operations is, under ideal 
conditions, a replica of the input FM signal accompanied by a relatively narrow 
band of noise surrounding the instantaneous frequency of the input signal. A 
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Fig. 7 Band-dividing scheme. 







conventional FM demodulator driven by this signal will have a lowered thresh- 
old because of the reduced noise power at its input. 

The amount of threshold reduction afforded by such a scheme depends 
upon two major factors: 

(1) The ability of the estimator to furnish an 
error-free estimate of the input frequency 
within a specific accuracy in the presence 
of noise, and 

(2) The presence of a coherent difference fre- 
quency signal which can be narrowband 
filtered without introducing serious fluctua- 
tions in amplitude. 

Each of these requirements must be carefully examined. In what 
follows, we take up the problem of estimation for the various possible forms in 
which the desired baseband signal can be presented to the FM demodulator at 
the transmitting end. The coherence requirement is now under study and will 
be covered in a later report. 

3 - 4. 1 Estimator Requirements for Various Message Representations 

The band-dividing concept involves a determination of the input 

instantaneous frequency by a set of contiguous bandpass filters or some other 

3 4 5 

equivalent structure. ’ 2 3 * 5 There are at least three ways of processing the 
message prior to frequency modulation: 

(1) No processing other than pre -emphasis by a 
linear filter, 

(2) Sampling and quantization into a set of fixed 
levels, and 

(3) Sampling and holding with no quantization. 
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The required number of filters and the filter bandwidths are different for 
the three cases. 

Consider first the case (No. 1) of an analog presentation 
of the baseband signal to the FM modulator. At the receiver, the "signal 
estimator" is made up of a suitably designed band contiguous bandpass filters 
each followed by an appropriate response measuring system. The important 
point here is that the (analog) instantaneous frequency must stay within the 
bandwidth of a given demodulator filter long enough for the filter output to 
build up to a response closely approximating the steady state, otherwise 
noise and distortion will dominate. If the filter bandwidth is given by B 

a, i 

then the instantaneous frequency must remain within B . for no less than 

I 

27r/B sec. Assuming a worst-case situation of sinusoidal FM of index 
3. f X 

6 at the maximum message frequency, co , the instantaneous frequency is 
given by 

w.(t) =co + 6 to sinco t. (3.47) 

l c m m 


The maximum instantaneous frequency rate is, 


dw. 

l 

dt 


= radians/sec^. 

m 


(3. 48) 


max 

Thus, the value that must be exceeded by the filter bandwidth is determined 
from 



max 


At 


6co 2 (^) < B 
m B — a,f 

9. y X 


(3. 49) 
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whence 

B a jf > co m /6 . (3.50) 

This simple result shows that the filter bandwidth in general must be con- 
siderably greater than the Nyquist rate of 2w and that the excess bandwidth 

m 

increases with the modulation index. A similar relationship holds for the 

optimized second order phase-locked loop, and, as was shown in Section 3. 1, 

for the usual i-f amplifier followed by a conventional FM discriminator. On 

this basis, the band-dividing approach clearly does not offer a worthwhile 

approach to FM demodulation. 

Consider next the case (No. 2) of sampling and quantizing 

prior to frequency modulation. At the receiver each level of quantization has 

a corresponding filter. The allowable minimum sampling frequency, f , is 

s 

twice the message bandwidth. Thus, if the highest significant frequency in 

the message spectrum is slightly less' than f cps, then the widest allowable 

time spacing between the periodically timed samples is 1/f = 1/2 f sec 

s,min m 

Consequently, the radiated FM signal is made up of a sequence of rectangular 
pulses 1/2 sec wide, each of which envelopes an r-f signal having one of a 
finite set of frequencies. It is well-known that such a frequency-hopping signal 
in the presence of white , gaussian noise is best estimated by means of a bank 
of filters each matched, to one of the rectangular r-f pulses. The normalized 
magnitude squared of the transmittance of each of these filters (shifted down 
too = 0) is given by 
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(3. 51) 


— AD / 


C (JSi 


[ 


sin (co / 2f ) 

-] 


CO 


/2f 


The (noise) bandwidth of such a filter is given by the integral of (3. 51) between 
-oo and +co, which yields 


B = 27rf rad/sec. (3. 52) 

m,f s 


For the allowable minimum value of f , 

s 


f . = 2f 

s,mm m 


the bandwidth B „ of each matched filter is 
m,f 

B . . = 2tt * 2f (3/53) 

m,f,mm m 

Since number of levels necessary to obtain a low level of quantization noise is 
usually large, a large number of filters spaced at least by twice the message 
bandwidth are required. This system therefore must have a deviation equal 
to the number of quantization levels times the message bandwidth, or a modu- 
lation index of N = number of levels. Usually the required index will be imprac- 
tically large, which makes case No. 2 not particularly interesting. 

The third case (No. 3) is demodulation of a signal in which the 
message has been subjected to a sample-and-hold operation prior to frequency 
modulation. Here, the instantaneous frequency is constant between samples 
and the filter bandwidth can be equal to the sampling frequency. The sample- 
and-hold operation does not destroy the analog message provided that the sam- 
pling rate is greater than twice the message bandwidth. 
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With no quantization employed, the constant values of fre- 
quency between pulses can assume any of the values in a continuum that covers 
the range of expected frequency excursions. Thus, the advantage of pre -fixing 
a b ank of filters matched to a discrete set of pulses — made possible by quantization 
— vanishes. Only sub-optimum solutions are now practicable. One such solution 
is to use abank of filters that are matched to a discrete set of the frequencies in the 
desired range. The specific selection of frequencies maybe made either on an ad 
hoc basis to be uniformly spaced within the desired range, or to be distributed on 
the basis of the probability density function of the baseband message process thus 
favoring the more probable frequency positions. One may also consider a dis- 
crete set of adaptive matched filters each of which acts like a phase -locked loop 
that can be frequency pulled only within a designated frequency sub-interval. 

The performance . of the estimator in the present case (No. 3) 
can be analyzed on the basis of determining the probability of error in identify- 
ing the frequency position of the signal from observation and comparison of the 
filter outputs. Such an analysis should take proper account of the fact that the 
filter most nearly matched to a particular signal pulse will be detuned relative 
to the actual pulse frequency. The amount of detuning will be distributed in 
accordance with the conditional probability distribution of the message sample 
values in the designated range. The effects of "matched" filter mismatch and 
detuning upon the probability of error are now being studied and will be covered 
in a future report. 
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Akima has considered only the demodulation of a signal 
frequency modulated by a sample -and -hold representation of the baseband 
message. His results suffer from two significant shortcomings of the analy 
sis: 


(a) . He erroneously applied error probability formulas, 

derived by Reiger° for the proper matched filter 
case, to the present case in which the filters can 
be detuned from the actual frequencies by up to one 
full baseband bandwidth. 

(b) He erroneously assumed that under high SNR con- 
ditions the filter bank will estimate the instantan- 
eous frequency within + &> where is the mes- 
sage band limit. This assumption ignores the 
potentially significant FM transient effects. 1 

Something must further be done to refine the frequency esti- 
mates after the estimator; otherwise, there will be a quantization noise which 
is unacceptably large for FM signals of reasonable modulation indices. Akima 
proposes to measure the instantaneous frequency at the output of the band- 
dividing filter having the largest amplitude. His handling of this process and 
assumptions about its results are inadequate. While this measurement, in 
principle, could be accomplished, the amount of circuitry involved would be 
rather large. A more workable alternative is the system proposed here in 
Fig. 7. 


Near the threshold, the estimator in Fig. 7 will occasionally 
make a mistake in selecting the proper filter because the noise amplitude in 
some filter will exceed the amplitude of the filter which contains the signal. 
Under moderately high SNR conditions, the probability of a mis -selection is 
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given by (N-l) times the probability of error for the corresponding binary 

g 

symmetric system (2 filters). These mis-selection errors constitute the 
threshold noise (or "error noise" as defined by Battail) of the band dividing 
demodulator. While the error rate ‘will depend upon the sampling rate and 
to a certain extent upon the number of channels, for a given sampling rate 
and about 10 channels the threshold noise will become significant when the' 

SNR in a bandwidth of 20)^ drops to about 12 db. (Akima’s intrinsic SNR is 
measured in a bandwidth of co^ and is therefore 3 db higher). 

Therefore, regarding the selection of a representation scheme 
for the input analog message it is evident that the sampling scheme which per- 
mits the smallest possible bandwidth for the estimating filters will give the 
lowest threshold. Also the sampling scheme must not introduce an excessive 
amount of quantization noise. Of the sampling schemes considered, the sam- 
ple -and -hold operation is best in this regard since it gives a minimum estima- 
tor bandwidth of 2co m regardless of the modulation index. Straightforward 
analog modulation requires an estimator bandwidth larger than this with the 
bandwidth increasing as the square root of the modulation index. 


41 

ADVANCED COMMUNICATIONS • RESEARCH AND DEVELOPMENT 




3. 5 Effects of Non-matched Filtering and the Interpulse Frequency 
Measurement in Band-Dividing Estimators 

Assume a sample-and-hold message frequency modulated on a 
carrier. The estimator contains a bank of filters, each filter having an 
impulse response described by, 

h (t) = 2g(t) cos a> t (3. 54) 

n 43 n 

where g(t) is real and lowpass . The filters are not necessarily matched 
to the input signal pulses. During any one signal pulse, the filter is hit 
with an input consisting of 

r(t) = s(t) + n(t) ; t > 0 

= n(t) ; t < 0 (3. 55) 

where s(t) is the pulse, of signal and n(t) is white gaussian noise. The output 

of the filter at the end of the square signal pulse is 

T T 

e (T) = r s(t) h (T-t) dt + f n(t) h (T-t) dt. (3, 56) 

o J n J n 

o -oo 

Adopting complex notation the output is 

T j(w -gj )t jto t 

e o (t) = j J g(T-t) e S n dt + V n 1 e n (3, 57) 

o 

where has the magnitude and phase of the (narrowband) filtered noise in 
(-3. 56), The envelope of the filter current is the magnitude of (3. 57).. The magnitude 
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of the signal termin(3.57)is essentially the magnitude of the Fourier Transform 
of the lowpass impulse response from 0 to T evaluated at co g - or 


l S out (T, l * 


g(t) 


-jut 


dt 


CO - Ci) 


'CO 


(3. 58) 


n 


There will be a degradation in (3,5 8) over the matched filter case if g(t) is not the 
impulse response of the matched filter and because of the detuning effect. The 
probability of error in selecting which filter has the signal will therefore in- 
crease. If the inter- channel interference problem is ignored, the degraded 
S/N ratio simply enters into the usual expressions for the system probability 
of error since the decision problem remains unchanged. The detuning is, of 
course, a function of the message and an exact statistical description of the 
detuning is difficult to obtain. Related to the detuning is the inter-channel 
interference caused by the response of more than one filter to the signal pulse. 
This effect will increase the probability that near the threshold a channel adja- 
cent to the signal one is erroneously selected. With a signal term present in 
more than one filter the decision problem itself is modified and available re- 
sults no longer apply. 

Results for Single -Pole Filters 

The matched filter is optimum in the sense that it maximizes the ratio 
of peak output signal power to average output noise power for white noise and 
a given signal waveform. For any filter, an efficiency factor may be defined 
as 
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_ peak output signal power 
a output noise power 


(3. 59) 


for a given ratio of input signal energy to noise power density. In the band- 

dividing demodulator, the input signal is a square pulse of width r and unit 

s 

amplitude. It is convenient to define the filter's zero-frequency response as 
2 

I H( jo) f = 1. The output noise power from the filter is given by 
-v N 


00 


n 


out . 2ir 


|H(jco)| dco 


-oo 


N 

if B n 


(3. 60) 


where is the usual 2-sided noise bandwidth of the filter in rad/sec-. Set- 


ting N o /2?r arbitrarily equal to 1, the efficiency factor becomes 

T 


n(t) dt 


n = 


/B 
V n 


with the additional restriction that 


(3. 61) 


! h(t) dt | = |H(jo)| = 1 . 


-CO 


The matched filter for a square pulse has a transfer function; 

sin «r / 2 

H <» = 


m 


a noise bandwidth of 


cor / 2 
s 


B 


2jt 


n, m 


(3. 62) 


(3. 63) 


(3. 64) 
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and a peak output signal amplitude of 1. Its efficiency factor is thus 


n 


m 


s 

2lt 


= 0 




C3. 65) 


The single-pole filter has an equivalent lowpass transfer function 


a 


H (jw) = — . 

0! J Jco + O. 


(3. 66) 


a noise bandwidth of 


B - i xa , 
n, a 


(a. 67) 


and an impulse response h^ (t) = a e 
Its efficiency factor is. 


-a t 


n 


T s 

1 - e s 


a. 


J net 


( 3 . 68 ) 


The efficiency is maximum when ct = 1. 25/T s , and is given by 


n 


a 


max 


= 0.36 f T~ . 
V s 


(3. 69) 


The single-pole filter is therefore worse than the matched filter by 
20 log (0. 4/0. 36) = 0. 9 db. 

The 3-db half -bandwidth of the most efficient single-pole filter is 


f = = 0. 20 f 

CO 2?r £ 


(3. 70) 


where f is the sample rate. It is obvious that even under steady-state 
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conditions a spacing of these filters by f (the sample rate) -will lead to 
"holes" between filters where the signal response is several db down. 

There are two solutions to this problem. One is to reduce the spacing 
between filters and add' more of them. The other is to increase the filter 
bandwidths, accepting some further sacrifice in efficiency. A factor of 
2 increase in bandwidth will reduce the efficiency by 1. 4 db. (The reduc- 
tion in efficiency is less than 3 db because the noise power doubles and 
the signal output also increases somewhat. ) The solution of adding more 
filters appears to offer less performance degradation. The probability of 
error in selecting a filter goes up roughly as the number of filters but a 
fraction of a db increase in input SNR will compensate for a 2 or 3 times 
increase in the number of filters. Barring practical considerations of com- 
plexity, it appears that a filter spacing of f /2 or f /3 is best for single- 

s s 

pole filters. The estimated frequency will then be quantized to finer steps 

and errors caused by inter-channel interference will be restricted largely 

to 1/2 or 1/3 off instead of f . 

s s 

It might be well to point out that even matched filters spaced by f 

s 

will be 4 db down at the overlap point, and closer spacing of these filters is 
also necessary. The spacing of f g suggested by Akima is based on the argu- 
ment that the bandwidth of a pulse r_ long is roughly ± f = 1 / t . The 

® s s 

question of the optimum filter spacing is unresolved and cannot be obtained 
from Akima' s model. 
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IV. EVALUATION OF CERTAIN PUBLISHED ANALYSES 


4. 1 Quasi -Linear Model of PLL Demodulator 

1 2 

In two recent papers J on phase -lock reception, Develet has advo- 
cated a "quasi -linear" model of phase-locked loops for obtaining "a simple 
analytical threshold criterion. 11 The proposed model and the treatment and 
results presented deserve careful assessment and criticism. 

We start with the first paper, which considers "reception of an 
unmodulated sinusoid embedded in additive white gaussian noise . " The 
discussion starts with a "simplified block diagram of general phase -lock 
receiver with no signal modulation imposed" in which the input noise is 
represented only by its quadrature component denoted Y(t). That this model 
and representation are questionable may be immediately seen from Eq. 

(3. 7 ) with ib( t) = 0, namely 

£ Lc (t) - {- l 1 + ^ X c, if (t) ] Sln hsc (t) 

+ B 3 - X q, if (t) COS h S c (t, } S hp (t) < 4 ' l > 

s 

subject to the condition |e(t)| = |<J> (t) | < w/2. The model suggested by 

OS c 

this equation is shown in Fig. 8. The corresponding model presented by 
Develet in his Fig. 2 differs from the correct one based on Eq. (4. 1) by the 
unjustifiable omission of the block marked cos ( ) in our own Fig. 8. 
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two multipliers within the loop, and the [^1 + ~~ x^ ^(t) J component of the 

s 1 

input. These omissions are objectionable because: 


a) The elimination of cos x imposes a more stringent 

restriction on |x| than does the linearization of sin x - 

the first rejected term in the approximation cos x 1 is 
2 3 

x /2 as opposed to x /6 for sin x «x. Therefore, omit- 
ting cos x while at the same time retaining sin x unlinear- 
ized is indefensible and not allowable. 

b) The omission of the input component 1 + — x .„{t) 

Hj g C,1I 

amounts to approximating this component by unity. This 

approximation is allowable provided that ~ |x (t)|«l 

ili g c,u 

almost all of the time and this is practically considered 
to hold only if 


(S/N) > 5 (or 7 db) 


With regard to (b), note that even if one goes out of the way and sets 

sin <[> (t) « <{> (t) and assumes that the closed-loop noise bandwidth of the 

osc osc 

resulting linear (not quasi -linear ) model is much smaller than the i-f band- 


width (which is usually justifiable for the scope of Ref. 1, but not in general 

for Ref. 2), the requirement for omitting 1 + x ,„(t) from the model is 

s c ' lf 

only reduced to 


(S/N) cjg > 5 (or 7 db) 


in which (S/N)^ is the ratio of mean-square signal to mean-square noise 
as filtered by the bandpass analog of the closed -loop low -pass equivalent. 
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Since the object of Develet's treatment is the analysis of a threshold which 
according to his computations should occur down at a S/N within the closed - 
loop bandwidth of 1.34 db, the indicated omission is completely unallowable. 

The quasi -linearization and analysis of the model considered by 
Develet in his Fig. 2 are also questionable. The quasi -linearization pro- 
cedure employed "essentially determines the average gain of the nonlinear 
device [sin ( )] under the expected operating conditions. " The objections 


against this are : 

a) The replacement of sin ( ) by an average gain is subject 
to restrictions on the allowable values of the argument 
of the sine function which are neither stated nor heeded 
by Develet. 

b) This replacement is definitely invalid for a phase -locked 

loop unless i (t) (or e(t) in Develet l s notation) is re- 
osc 

stricted in magnitude, because it cannot account for the 
break of lock that results whenever the condition 

|<j> (t) | ^ v 12 

is violated. The condition for lock can be expected to be 
satisfied almost all of the time as long as 

a 2 = <j> 2 (t) < tt 2 /6 4 « 1 /6.4 

where we have allowed a gaussian-like character for 

<j) (t) with a crest factor of 4. This shows that Develet's 

osc 

model breaks down for a > 7r/8 0. 4 rad. Account of this 

breakdown of the model must be taken in defining the limits 
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of the integration in Develet's Eq. (2). But since the 
results and curves cannot be used for a > 0. 4 rad, the 
necessary correction for the limits of the integration 
will not have a substantial effect. According to the 
curves of Fig. 5, the "loop signal to noise ratio" for 
q- = 0. 4 is 2 db for the "quasi -linear approximation" 
and 5 db for the "completely linear approximation. " 

These are then the thresholds that can be read from 
Develet's computations. (A more complete account 
of this threshold is presented in Ref. 3. ) 

Philosophically, the concept of replacing "instantaneous" effects by 
the effect on the average is disquieting in the present application because 
the breakdown of performance is caused by disruptions resulting from ex- 
cesses by the instantaneous values of the phase error that cause the system 
to be either in lock , and hence describable by the model used, or out of 
lock , and hence not describable by the model used. 

It also appears that approximation of a substantially nonlinear 
input-output characteristic by its average gain over the distribution of an 
input random process is a great oversimplification in the computation of 
the effects of this nonlinearity upon the input process, or upon the sum of 
this and another process. 

It is of interest to point out that the curve marked "quasi -linear 
approximation" in Develet's Fig. 4 is incorrect. The reason for this is 
that the in the left-hand member of Develet's Eq. (7) is the noise 
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bandwidth of the closed-loop system of his Pig. 3. This B is therefore 

n 

function of or. In plotting his curve, Develet has ignored this dependence 
on the basis of the statement that holding B n constant ". . . can theoretically 

9 

be accomplished by varying the gain of F(s) in opposition to exp (-a /2). " 

2 

Varying the gain of F(s) in opposition to exp (-a /2) means precisely replac- 

2 

ing F(s) by F{s) exp (a /2). If this is indeed done, the open-loop transfer 
function of Develet 1 s model becomes identical with that of what he calls 
the "completely linear" model, and the "quasi -linear" model therefore be- 
comes "the completely linear" model. Ignoring the dependence of B upon 
o is therefore unjustified and the "quasi -linear approximation" curve in 
Fig. 4 does not correctly portray the performance of the model of Fig. 3. 

To sum up the above critique of Ref. 1: 

a) The model of Fig. 2 is an incorrect representation of 
the "general phase -lock receiver" described in Fig. 1. 

b) The model of Fig. 2 cannot be used for analysing the 
performance of phase-lock receivers at or near their 
threshold of performance breakdown. 

c) The quasi -linear model employed omits vital effects 
of nonlinear behavior at or near the threshold of per- 
formance breakdown. 

d) Even when proper corrections are made throughout for 
model, expressions and curves, the validity of the 
results and curves must be restricted to values of rms 
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loop error a less than approximately 0. 4 rad, because 
above this value of a a quasi-gaussian phase error 
distribution would indicate that unlock becomes likely 
during an intolerable fraction of the time. 

2 

In his second paper, entitled, "A Threshold Criterion for Phase - 
Lock Demodulation, " Develet considers phase or frequency modulation by 
"Gaussian signals corrupted by additive white Gaussian noise. " At the out- 
set, a restrictive assumption is made and used in modeling the system 
which restricts the validity of the model only to the case in which the pre- 
detection bandwidth is much larger than the phase -lock loop bandwidth. This 
assumption makes it allowable to consider the noise in the low-frequency 
output of the product demodulator to have a gaussian character and a cor- 
relation function (or spectral distribution) identical with that of the cophasal 
and quadrature components of the input noise. (The results of the analysis 
are applied extensively in a later portion of the paper to cases in which this 
basic assumption does not hold at all. Specifically, note that curves are 
plotted and discussed for rms deviation ratios of less than 10 in Figs. 4 
and 5. ) As in the first paper, discussed above, the fact that break of lock 
ultimately limits the usefulness of the output of a phase -locked demodulator, 
and hence incurs the ultimate threshold, is completely ignored in this 
second paper. Thus, the basic model and analysis used in this paper are 
inadequate and invalid because: 
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a) The error signal e(t) is not restricted between ±tr/2 
in order to avoid unlock. 

b) Unlock introduces severe disturbances that cannot 
be characterized by the model of. injected noise used 
in the paper. 

c) The replacement of the nonlinear characteristic, sin ( ), 
by its average derivative oyer a gaussian distribution is 
improper because the error signal e(t) is not in general 
gaussian, and because important by-products of the non- 
linearity are thereby excluded whose instantaneous 
(rather than average) fluctuations have a critical effect 
upon the threshold performance of a phase -locked loop. 

d) It is not realized in the formulation of the model that 
in order to keep nonlinear distortion in the absence of 
noise below tolerable limits, the magnitude of the signal 
component of e(t) must be restricted to values well be- 
low unity. 

4.2 D eyelet's Least Squares PLL Optimization 

2 

Develet has also carried out a least squares optimization analysis 
of the PLL demodulator. There is doubt about Develet 1 s results for two 
reasons : 

2 

1. He uses an unrealistically large loop error, o =1. 

2. He uses the quasi-linear equivalent gain model which 
seems to be of little value because it considers only 
one of the effects of the nonlinearity. 
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These doubts justify repeating the analysis with a considerably smaller 
2 

value of cr (0. 1) and without quasi -linearization. However, the results 
of our analysis (Chapter V of this report) agree with Develet on the one 
point most important to this project: there is about a 3 db difference in 
threshold between an optimized second-order system and the optimum 
system at a modulation index of 1 0. 

4. 3 Sanneman and Rowbotham Paper 

The digital simulation of a second-order PLL by Sanneman and 
4 

Rowbotham has produced an interesting result which was not fully devel- 
oped in their paper. They used a sampled -data system to simulate the 
PLL; the sampling rate being chosen large enough to ensure correspond- 
ence with the continuous system. The computerized loop was supplied with 
a simulated carrier plus noise input and allowed to run (starting with pre- 
scribed initial conditions) until a break of lock occurred. A large number 
of runs were taken for each set of conditions so that an average time to 
unlock could be computed for various carrier -to -noise ratios and initial 
conditions. 

The authors make an issue of the fact that the probability of unlock- 
ing in T seconds can be described by the empirical law 

-co f/oo T 

P(uy) = 1 - e n n , (4.2) 

where T is the average time to unlock. This function is simply the probability 


of one or more Poisson events in an interval to j and would be obtained for 

n 
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a model where the probability of -unlock is constant in any given tim e inter- 
val. In other words, their data supports the hypothesis that unlock events 
occur individually and collectively at random. 

Of greater interest is the plot of average time to unlock versus 
phase error given in their Fig. 12. The abcissa of this plot is directly the 

rms phase error and the ordinate is the log of the normalized time to un- 

3 

lock. Viterbi has obtained an expression for the mean time to unlock for 
the filter -less PLL which takes the form 

— 7 r f 2. 0 “1 

B r « i exp [— ] (4. 3) 

where is the noise bandwidth of the loop. A plot of log B f versus 
2 

1/^ yields a straight line. This result immediately suggests a re -plotting 
of S and It's Fig. 12 with an abcissa of 1/ . Such a plot is shown in our 
Fig. 9. The data now yields a straight line. For both initial conditions 
equal to zero, S and R*s data is described by 

= 2exp[^~-] (4.4) 

Equation (4. 4) therefore appears to be an adequate empirical expression 
for the mean time to unlock in the second -order PLL, a result which Viterbi 
was unable to obtain theoretically. The use of Eq. (4. 4) is further justified 
by the fact that it is very similar to the theoretical result in Eq. (4. 3) for 
the filterless loop. 
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Fig. 9 Replot of data in Sanneman and Rowbotham's Fig. 12. 
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It is now possible to define a threshold value for in terms of a 

prescribed number of unlocks per second. For example, a second -order 

loop having a 16 kc natural frequency (co^ = 105) will have an average of 1 

2 2 

unlock per second at a mean-squared phase error of cr^ = 0. 15 rad . For 
2 12 

= 1 as proposed by Develet J and others the corresponding unlock rate 

is 10,000 per second for co^/2?r =16 kc. 

These results are for an unmodulated carrier plus noise and zero 

steady-state phase error. If the mean-squared phase error includes a 

random modulation-tracking error the unlock rate will certainly be no less 

than that given by Eq. (4. 4). Therefore, for audio -type loop bandwidths 

and a threshold defined by an unlock rate of about 1. second the threshold 

2 

occurs at a mean-square phase error of 0.15 rad . The unlock rate is 
much more sensitive to the phase error than it is to the loop bandwidth, 
therefore one need not be extremely precise in specifying the loop band- 
width or the exact unlock rate. For example, a 10 times decrease in unlock 

2 

rate corresponds to = 0.123 instead of 0.15 or about 0.9 db difference 

in SNR measured in the loop bandwidth. 

On the basis of Sanneman and Rowbotham's results, it therefore 

appears that an appropriate threshold value of total mean -squared phase 

2 

error in the second-order PLL is 0.15 rad . At this level of error, the 
linearized model of the PLL is reasonably accurate. The average gain of 
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the multiplier, as computed by D eyelet"*", is 0. 93 instead of 1. The threshold 


criteria of Develet 5 and Van Trees which are based upon a point of 
divergence of me an -square phase error in an approximate model yield un- 
reasonably high rates of unlock. The divergence of phase error (or gross 
unlock in some sense) is a phenomenon associated with the particular model 
being used. This phenomenon of the model is not an adequate description 
of the random breaks of lock which are observed to occur near the threshold. 
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V. HIGHER -ORDER PHASE -LOCKED LOOP PERFORMANCE 
5. 1 Introduction 

The use of a sinusoidal FM message for the purposeof comparing 
different PLL structures is supported by the following reasons: 

a) The analytical evaluation of system performance is rather 
straightforward which permits one to keep a close tag on 
the analysis assumptions and limitations, common causes 
of serious trouble in the study of PLL systems. 

b) A comparison with existing conventional demodulator anal- 
yses is possible., though the PLL threshold analysis must 
demand both message and noise to be present at the same 
time due to the signal tracking vs. noise rejection compro- 
mise of the system, a condition which may be somewhat 
relaxed in conventional demodulation. 

c) The demand to match the laboratory experiments to the 
theoretical analysis can be met more easily. 

Meanwhile, the use of higher order loops is perhaps motivated by 
the work of Jaffe and Rechtin* where a second order loop was found to be 
optimum for a frequency step message and a third order loop for a frequency 
ramp message. They considered a linear, noisy (additive noise with flat 
power spectral density), time invariant model and a minimum total (message 


* R. Jaffe and E. Rechtin, "Design and Performance of Phase-Lock Cir- 
cuits Capable of Near -Optimum Performance over a Wide Range of 
Input Signal and Noise Levels", IRE Trans, on Info. Theory, March, 1955. 
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distortion plus noise) mean-square error as the optimization criteria. In 
any case, the only hint we want to draw is that the third order loop can be 
expected to do a good job on the small angle region of the sinusoidal FM 
signal where sinx & x allows a ramp approximation, while there is no 
equivalent motivation for a frequency step characterization and, conse- 
quently, a second order loop. 

5. 2 Phase Error and Threshold Considerations 


The linear, time invariant phase model of a PLL is assumed as a 
starting point. This inherently assumes an approximation to the problem, 
since the loop phase noise and the additive input noise will show some 
statistical peaks that exceed the linearity and time invariancy- conditions 
with a finite probability such that the model in question is not valid always. 
However, such behavior is assumed to have a small frequency of occur- 
rence so that the model may be assumed as valid "almost all the time. " 
The threshold phenomena essentially represent the breakdown of this 
model and are thus characterized by the occurrence of the statistical peaks 
with a frequency (probability measure) larger than allowable. 

Consider the model of interest as shown in Fig. 11, where the 
input phase is the message phase cj)^ plus the normalized quadrature 
noise <j> , and 4 > vco ^ K and F(p) respectively represent the YCO output 
phase, loop gain and loop filter. The loop phase error in tracking the 
message <j> is given by 
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K {s) = ( f ) s (s) - Kco {s) = ^S (s> ■ H(s)[< j ) s (s > + ( f > n (s) ^ 

= [1 - H(s)]<j> (s) - H(s)cj> (s) (5. 1) 

s n. 


where the first term is a message distortion term characterizing the 

inability of the loop to exactly reproduce the message and being denoted as 

the (linear) distortion error tj) , , while the second term is a noise term 

e,d 

representing the VCO phase noise caused by the additive input noise and 
being denoted as the noise error <j) 

0 J 13 

In the case of a sinusoidal FM message having a modulation index 6 
and a modulation frequency « rad/ sec, i. e. , 

<j> (t) = 6 sinco t , (5. 2) 

T s m 

the distortion error is sinusoidal with an amplitude given by 

^e,dl = 6 1 X “ H Cj“ m >l rad (5>3) 

so that j 1 - H(jco^)| essentially represents a compression factor for the 
modulation index. Meanwhile, if the noise has a normalized flat power 
spectral density of <3?/2S (rad/cps) and zero mean, the rms phase noise 
error is given by 

i = rJL f°° i m - u 2 ^ - ] 1 / 2 

^e,n (rms ) L 2S J 2t r J 

<|>B v. 1 / 2 

= (rdr) rad (s - 4) 

oo 2 

where B = J j H( j 2wf) | df is the equivalent phase noise bandwidth of the 

-CO 
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system in cps. If the statistics of this phase noise are known, then the 
probability of exceeding F times the rms value can be evaluated for any F 
and a statistical peak is thus defined when a sufficiently small probability 
is achieved. 

On this basis, the total peak phase error may be taken as 



I* 


e,d 


+ F<b 

e ' n (rms) 


(5.5) 


where the "crest facto'r" F may be readapted further to acknowledge the 

fact that the distortion and noise error peaks do not occur necessarily at 

the same time instant. Finally, if <j> characterizes the critical phase 

c 

error beyond which the model is assumed to break down, then operation 

above the threshold "almost all the time" is characterized by I < <i> . 

The choice of A is linked to the crest factor F. If the noise 
T c 

statistics are assumed gaussian, then the use of F = 1, 2, 3 implies that 

the "peak" value F<|> will be exceeded only 33%, 3%, 1% o'f the 

e,n(rms) 

time, respectively. On this basis, the use of F = 3 and <j> c = 7r/2 seems 
suitable. In reality, the noise is not gaussian near threshold and the PKL 
model is nonlinear near jt/ 2 and time variant near threshold. Still, these 
parameters are but a performance criterion and should provide an adequate 
guidance to the threshold phenomenon. A proof of this statement is perhaps 
the fact that 3^^ < tt/ 2 (neglect distortion temporarily) implies <j> < 

tt/6 which is a common upper bound in the linear analysis of PLL systems. 
The inclusion of the distortion term will yield a somewhat smaller rms 
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noise error but not by a considerable amount if the loop is operating in an 
optimized behavior , as shown in the sections that follow. 

5. 3 The Second Order Loop 

The second order loop is characterized by a loop filter of the form 
1 + TfS 


F(s> = 


1 + r o s 


<5. 6) 


which results in a noise bandwidth given by 


KTi 

1 + — - 
K T 2 

B n " 2 »i;:, ° PS 




1 /• 1 Kt 1-\ 

~Gt + tv) for KT 1 >:> 1 


(5. 7) 


At this point it is useful to notice that this bandwidth uniquely 
specifics the noise error for a given input noise density <1>/S (and noise 
statistics) so that several second order loops having different design 
parameters K, * r ^ but a similar will show the same noise rejec- 
tion capability and their relative merits will only be determined by their 
message tracking (distortion error) capabilities. On this basis j it seems 
useful to change the set of design parameters (K, t to a new set 
that has as one of the parameters. A look at the approximate relation 
in Eq. (5. 7) suggests the following transformation: 

1 KT 1 

— = aB , — ~ = (2 - a)B , 0 < a < 2 (5.8) 

T X n r 2 - n 
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which yields the new set of parameters (K , a) where* 


aB 


t 2 = 


n 


K K 

a(2 - a) B * n 
n 


» a 


(5.9) 


The peak distortion error in terms of the new parameters is given, 
for the second order loop, by 


1 + 


e.d 1 


CO 


(if) 


+ 


a 2 (2 - a) 2 


CO 4 

(if) 


1 + 


2 - 3a 


co t. 

(if) 


a (2 - a) n 


a 2 (2 - a) 


co 4 

2 vB ) 
n 


1/2 


(5. 10) 


K. 


for r— 1 y> 2a, and by noting that the maximum value of a(2 - a) is 1 and 
B n 

that K can be made very large, the first term in the root numerator may 
be dropped. Also, it seems reasonable to expect that either the first or 
the third term in the root denominator will predominate, except perhaps 
for a small range of co m /B n values, which results in a compression factor 


2_V/2 


of the form ( ^ ' x ~ + x /) 
is to occur, i. e. , 


and x « 1 is desired if any noticeable compression 




6 co 


m 


e,d 


(5. 11) 


a(2 - a) B 


n 


* The exact noise bandwidth is actually given by B^/[l + (aB^/K)] 

approximates for K/B^ 2S> a, a reasonable relation since 0< a< 2. 
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Thus, optimum performance occurs for a = 1 where the compression is 
maximum for a given noise error (i. e. , B^). 


A more exact analysis can be done by plotting the compression 
factor with only the first term in the root numerator being neglected, i, e, , 

[<|> e d | « 6 [l + < 2 ~ a ><2 . l .3a) + a _. (2 _^a )— (5 . 12) 

e ' 7 y 


where y = w / B and the results are essentially shown in Fig. 11. For 
a constant y , i. e. , a given noise error, the interest is to operate at the 
largest possible abcissa in order to compress the modulation index as 
much as possible. Operation at y > 1 is inadmissible since no essential 
compression occurs, while operation at very small y 1 s is excellent though 
the price is an extremely large noise bandwidth and noise error. In gen- 
eral, it is noted that a = 1 is optimum for y < 1 where a compression 
occurs. Moreover, for y < 0.7 the approximation of the abcissa by y - 2 
(at a = 1) can be found to show a negligible error. On this basis, our 
choice of Eq. (5,11) with a = 1 is truly an optimum behavior. Notice 
finally that the neglected factor containing the w m /K term will not alter 
the a choice but w ill only introduce a reduction in the assumed com- 
pression. 

The total error after the first round of optimization (K large, a = 1) 
is thus given by 


6 CO 


4 > 


m 


e,p 


+ F 


.$B 


B 


V 2S J 


1/2 


(5. 13) 


n 
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and the final optimization is the choice of B^ for a minimum error, which 


can be found to be 


B = ( 

n opt ^ 


4 6u) 


2 2 / 5 


m-\ ,^2S\ 

) C -§) cps 


1/5 


F 


(5. 14) 


Moreover, if this error is to be kept below <|> c "almost all the time, " 
the threshold condition reads 


S . 1.75 F P 1/2 

$ * Ts J2 5 “m 




(5. 15) 


It is of interest to further analyze Eq. <5. 13) with given by 
Eq. (5. 14). Direct substitution will prove the fact that the noise error 
term will be exactly 4 times larger than the distortion error term under 
this optimized behavior. This means that the total error can be considered 
as 5 times the distortion error (the other 4 contributions representing the 
bandwidth limitation introduced by noise) or as 5/4 times the noise error 
(the other 1/4 contribution representing the distortion limitation). By 
using the first characterization, the threshold condition reads 


5 6"\ V 2 


! V * © 


w m °P S 


(5. 16) 


while the second characterization yields a threshold condition of 


(SNR) >0 

B “opt 


(5. 17) 
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These two conditions respectively give the loop bandwidth required for 
tracking the message in the presence of noise and the loop SNR required 
for maintaining message lock in the presence of noise, i. e. , the PL.L 
threshold SNR. 

If the intrinsic SNR is defined as 


(SNR). = <SNR) 0 . 

1 2f m 

then the intrinsic SNR at threshold is given by 
(SNR). = (SNR) 

From Eqs. (5. 16) and (5. 17) this ratio is 


n opt 

'B- ’ 2f 

n opt ni 


db (SNR). - 12 db + 1Q1 ° g6 

X £ 


(5. 18) 


for a crest factor F - 3 and <J> = ?r/2 . 

5. 4 The Third Order Loop 

The third order loop is characterized’ by a loop filter of the form 

2 

1 + T. s + (t„s) 

F(s) = 2 ~ (5.19) 

1 +r 2 s + (r s) 

which results in a noise bandwidth given by 


73 

ADVANCED COMMUNICATIONS • RESEARCH AND DEVELOPMENT 


1 + Kr x 




1 - 


K Va 


1 - 


Kr i T 3 




Kt : » 1 
Kt 3 2 >> T 2 


(5. 20 


The choice of a variable transformation in order to change the set 

(K, t j T , T , r ) to another one having B as a parameter is not 
i ^ o n 

obvious. However; an analogy with the second order loop suggests 


1 2 = 1 __ 

1 aB J T 3 , _ 2 
n b B 

n 



(5. 21] 


where 


JK 

B 


» a 


K 


n 


B 


» br 2 ; (a; b, c) > 0 


n 


2 ab c , b ,, 

and a + 2 - + — + — (5. 22) 

c b a 

On this basis, the new set of parameters is formed by K, B , T 0 and 
any two of the three parameters (a, b, c), the other one being defined 
by Eq. (5. 22). 


74 


ADVANCED communications 


RESEARCH AND DEVELOPMENT- 


AD, 


COM 


The amplitude of the distortion error is now given, for the third 
order loop, by 

r ^ 1/2 


IW 


6< 




i + -5-r— - r— - -^ + 2 

X ~ ' - 4 V 2 




k. 


2 k.,2 ~ a J 
y b 


b J 6 

y ^ 


, K _ K 

for — — » a , — 

n B 


2 >> bT 2 


n 


(5. 23) 


where y = W m /B n as before. Again it seems that K and r can be chosen 

to neglect the last two terms in the root numerator, and that the last term 

in the root denominator may be found to predominate for an effective 

compression of the modulation index. In this case, Eq. (5. 23) reads 

3 


I* 


6o> 


m 


^ cB 3 


(5. 24) 


n 


and the interest is to select the largest possible value of c such that both 
Eq. (5. 22) and the approximations that yielded Eq. (5. 24) are satisfied. 

In general, it seems that c > 1 can be obtained, though an upper bound 
is difficult to establish since it implies component feasibility considerations. 
On this basis, the total error is then given by 


♦ 


6 Cl> 


m 


4>B 1 / 2 


®' P cB 3 
n 


+ 


F C ^ 

V 2S ) 


(5. 25) 
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and the noise bandwidth optimization gives 


B 

11 4 . 

opt 


c 


6 6 co 

m 

cF 




cps 


<5. 26) 


and the threshold condition of the phase error being kept below <(> "almost 

T c 

all the time" reads 


S 

$ 


> 


1.29 F 

c 1 / 3 * 7 / 3 

T c 


.1/3 

0 co 


m 


(5. 27) 


An analysis of the total error analogous to that effected for the 
second order loop results in the noise error term being 6 times larger 
than the distortion error term in the optimized behavior. By considering 
the total error as 7 times the distortion error, the threshold condition 
yields the following tracking requirement: 

s 76 nI/3 

B > ( — t ) u cps (5.28) 

n . vcd) y m 
opt T c 


while by considering the total error as 7/6 times the noise error, the 
following loop threshold SNR is obtained: 


(snr) b 

n opt 



(5. 29) 


Corresponding to the second order loop the third order loop has an 
intrinsic SNR at threshold of 
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db (SNR). = <SNR) 2f = 11.1 db + 101 3 Og6 (5.30) 


for F = 3 , <|> c = 7r/2 and c = 1 . If larger values of c can. be obtained the 
SNR will be reduced by 101ogc/3 db. The threshold improvement over 
the second order loop is 0,9 + 1.7 log 6 db. 

5. 5 Systems Comparison 

A numerical comparison of the two system thresholds is dependent 
on the choice of the critical error <(> and the parameter c , as well as on 
the actual value of the modulation index. However, for <j> < w/2 arid 
c > 1 , the third order loop will always show a threshold improvement for 
reasonably large indices greater than one, the improvement being increased 
as c becomes larger. 

There is another important advantage of the third order loop. A 
look at Eqs. (5. 16) and (5. 28) shows that reasonably large bandwidths may 
result to such extent that the fact that they may become comparable to the 
i-f bandwidths must be accounted for. The third order loop is shown to 
have a smaller loop bandwidth for a given message so that its By,/B ratio 
will be larger. 

Finally, it is of interest to notice that the two loop thresholds, 
i. e. , Eqs. (5. 17) and (5. 29) are essentially similar. 
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5. 6 Higher Order Loops 

The results of the previous sections can he tentatively extended to 

higher order loops, i. e. , higher order loop filters. In general, if the 

k. 

distortion error term can be written as inversely proportional to , 

1/2 

the noise error being always directly proportional to ' , then the 

following relation will hold after the total error is minimized with respect 


to B : 
n 


2k + 1 

(noise error) = — — (distortion error) 
and the relevant effect will be the B„ and S/$ dependence on 6 


2k + 1 


such that both the system threshold and the required loop bandwidth are 
considerably reduced as k increases. The crucial issue is the capability 
of being able to express the distortion error in the desired form: on the 
positive side there are more parameters to play with while on the nega- 
tive side there are more approximations involved (the analysis is already 


difficult for the third order loop). 
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5. 7 


Least-Square Optimization of PL.L. Performance 


We shall next present the optimum performance attainable with a 
phase -locked loop and discuss implementation intended to achieve this 
performance. 

5. 7. 1 Goals of Optimization Procedure 

Tne PLL receiver of interest is shown in Pig. 12. The received 
signal r(t) is composed of the transmitted signal s(t) and noise n(t) which 
is white in the spectral region occupied by s(t). We make the following 
definitions: 

s(t) = *f~2 A sin(a> c t + 6m(t)) (5.31) 

m(t) = Message, a random process with 
spectrum S m (cj) 

6 = modulation index 

v (t) = '/2~ A cos + ^ *n^(t)) (5- 32) 

We decompose the noise as follows: 

n(t) = J~2 A n (t) sin (ut + S m(t)) + >J~2 A n (t) cos (w t + 6m(t)) (5. 33) 
1 C Ci c 

n (t) and n (t) are low-pass equivalents of n(t) and have a flat spectrum of 

X & 

2 

height N q /2A in the low-pass spectral region occupied by m(t). With these 
definitions we have (discarding double frequency terms which will not pass 
through the loop filter): 

r(t) x v vc0 (t) = A 2 Cl + n^t)] sin 6 [m(t) - m (t)] 

+ A 2 n 2 (t) cos 6 [ m(t) - m (t)] (5. 34) 
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Fig. 12 Model of PLL. 
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Thus we may draw the non-linear., time -variant phase transfer model of 


Fig. 13. There are two effects of importance here. First, if the loop error 
e^(t) becomes larger than tt/ 2 the loop will lose lock. Second, for small 
values of e^(t) the system is linear while for large values of e^(t) it is non- 
linear and time-variant. We can insure that operation will be linear most 
of the time by requiring that VAR[ e^(t)] < ty^ where is a suitably 
chosen constraint level. The magnitude -of VAR[ e (t) ] is also a measure 
of how often e^,(t) will exceed 7r/2, causing the loop to lose lock. The phen- 
omenon of threshold occurs when the loop loses lock with sufficient frequency 
to degrade performance. Threshold has no direct connection with the transi- 
tion from linearity to non-linearity, but will occur for about the same values 
of VAR[ e^(t)] which define this transition. Therefore, in the work that 
follows, we will postulate that threshold occurs for some loop error 
VAR[ e (t)] = a and that operation of the loop is linear until threshold is 

jL C 

reached. Thus, all following discussion will be based on the linear, time- 
invariant model of Fig. 14. Our goal will be to design F(s) so as to minimize 
VAR[ e^(t)l and thus optimize the system with respect to threshold. 

5. 7, 2 Linear Filter Optimization; Realizability 

It is important to recognize that we must require F(s) to be realizable. 
In an ordinary filtering problem we may approximate the performance of an 
unrealizable filter by introducing delay between the input and output. However, 
in the phase-locked loop the filter output is required in order to form its output, 
so we cannot use this stratagem. 


81 

ADVANCED COMMUNICATIONS • RESEARCH AND DEVELOPMENT 







ADVANCED COMMUNICATIONS • RESEARCH AND DEVELOPMENT 



AD, 


COM 


The problem of the optimum linear filter has been solved in a con- 
venient form by Yovits and Jackson^. Their results give the following 
expressions (in the white noise case) for the minimum error and for the 
filter which produces it: 


VAR[ V tn o Pt 


N .<» 
o 


2ttA 


2 J 


log 


1 + 2 -~ 6 2 S («)] dco 

N 0 m J 


( 5 . 35 ) 


and 


l-H(jto) 


1 + 


2A 

N 


6 S («) 
m 


( 5 . 36 ) 


where we have defined the closed loop transfer function H(jcu) by 


H(s) 


KF(s) 
s + KF(s) 


( 5 . 37 ) 


In order to calculate the loop error for non-optimum loop filters we must 

separate the two different causes of error. First there is a noise error e^ 

which is due to the presence of the noise source ^(t) 1°°P* Second there 

is a distortion error e, . due to the fact that the loop transfer function is not 

£d 

identically unity over the spectrum of the message. These errors are given 


by 


VAR[ e n (t)l = - 
Jin 7r 


® N 


§ | H< jto) | 2 dw 


o 2A 


( 5 . 38 ) 


and 


1 P 00 2 2 

VAR [ e (t)] = - I 6 S («) | l-H(jm) J dw 
*- led 7T o m 


( 5 . 39 ) 
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Because the message and noise are independent, these errors are inde- 
pendent., and we may write 


VAR[e (t)] = VAR [e (t)] +VAR[e <t)] 


5. 7. 3 Performance of Optimum Realizable Loop for PM with Rectangular 
Message Spectrum* ~ ~~ 


We consider the message spectrum given by 


co | < a. 


(5.41) 


co > O' 


We define the message bandwidth to be the bandwidth of a rectangular low-pass 
spectrum which has the same height at the origin as S (co) and gives the same 
total power. Thus 


If 

B “ o 77T\ | S ( co) dco 

eq S m (°) *Loo m 


(5.42) 


We will generally normalize S- (co) so that 

m 


7T~ i S (co) dco = 1 
07 r J m 


(5. 43) 


(Incidentally, this is the justification for calling 8 the modulation index), hence 


B - 2tt/S (co) 
eq ' m 

In the present case, of course, this reduces to 


(5. 44) 


B = 2a 
eq 


(5. 45) 


This is the approach taken by Develet, Ref. 2. 
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Therefore the signal -to -noise ratio in the message bandwidth is given by 


A = 


2ttA 

ffiN, 


(5. 46) 


Inserting (5.41 ) and (5.46) into Eqs.(5. 35) and (5. 36) we find 


VARCe /‘ )] opt ■ X (1-A6 2 ) 


(5.47) 


1 “ H(jto) i 


1 + A6 


u| < a 


(5,48) 


From Eq. (5.47) we can relate modulation index., message bandwidth SNR, 
and threshold as shown in Fig. 15. From Eq. (5.48) it is not clear how to 
construct the optimum loop filter; approximations to it will be discussed in 
section 5. 7. 5. 

5.7.4 Optimum Realizable Loop for FM with R-C Message Spectrum 
We take ' 


S m (w) " 2, 2 2. 

oj (to +a_ ) 


(5.49) 


and define 


(5.50) 


This spectrum is not normalized in accordance with Eq. (5. 43) because its 
integral is infinite. Therefore 6 is not properly a modulation index. 
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Then from (5. 36) we have 


1 ~ H( jco) j 


2 , 2 2 

CO (to +0! ) 


„ , s 4 2 2, *2 
1 + y 6 S (to) to ■+■ cl to +y o 
a 


Note that 


4 22 2 ,2 w 2 , v 

<o + a to +y 6 = {co + p ^ ) (to + pg) 


where 


P 1 = \ + y« 4 - 4v 6 2 ; P 2 = J (a 2 - J ol 4 - 4y 6 


So, 


1 - H(jco) j 


(— jco) (jco) ( g -t-jco) (Of - .jco) 


, 1/2 1/2 . 1/2 w 1/2 . , 

(p^ +]to)(p i - jco) (p 2 +Jco)(p 2 -Jco) 


and 


1 - H(jco) 


.jco(Qf + jco) 


, 1/2,. w 1/2 , 

(p i + 3co) (p 2 + jco) 


Whence, 


H(jto) 


(P x P 2 > lyr 2 + j“ (p ^ 2 + P 2 ly ^ 2 - » ) 
(P, 1 ^ jw) + 


Note p p = y 6 and furthermore 

JL u 


1/2 1/2 /, 1/2 1 / 2 / 
Pi + P 2 = v(P x + P 2 > 


/ 1/2 

Pl+P 2 + 2Pi ; 


/ o 2 1/2 

y Of + 2(y 6 ) 


(5. 51) 

(5. 52) 
) (5.53) 


(5. 54) 
(5. 55) 


1/2 

'2 

(5. 56) 
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So 


H(jw) - 


2 - 1/2 • ■ -L 

v O' 


+ jco U 


2 c 2> 

+ 2(y 6 ) 


1/2 


O' ) 


(p , 1/2 +W (p , 1/2 + j u ) 


C5. 57) 


The importance of this result is that it shows that the optimum loop for an 
FM system with an R-C message spectrum is a second order loop. That 
this is not true for PM' with a rectangular message spectrum is evident 
from the form of Eq. (5. 48). 


5. 7. 5 Approximations to the Optimum Realizable Loop Filter for a 
Rectangular Message Spectrum 

We found in section 5. 7. 3 that the form of the optimum loop filter is 
given by 


I 1 " H(jco) 1 2 


r 1 
1 + A 6 

v 1 


GJ < Oi 


o>j > a 


(5. 58) 


It is clear that to build a system having exactly this transfer function would 
require an infinite number of poles and zeros. However, we can certainly 
approximate it as closely as desired by using a system of high enough order. 
Because of the increase of complexity associated with higher order filters 
we would like to keep the order, n, as small as possible while meeting a 
desired performance criterion. One way of approaching this problem might 
be to find the "best" n^ 1 order filter by the following procedure. First, we 
would write the unknown filter F(s) in terms of 2(n~l) parameters as 
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n-l S + a . 

F(s) ■ TT srr: 

i= 1 i 


It Is easily verified by Eq. (5. 3 7) that this gives an H(s) with denominator 

of order n. Next we could vary the 2(n-l) parameters a, o'. -.a: , ft , 

i* n~JL l 

ft 2 ’ • • • & n -\ so as minimize the resulting loop error variance as computed 
by Eqs. (5.38) through (5.40). This procedure would, of course, be impossible 
to carry out for n reasonably large. An alternative approach is to assume 
more of the structure initially, leaving less parameters to be varied. 

Adopting the latter procedure, we notice that [ 1 - H(jco) | is equal to 
the difference between a constant and a rectangle: 


l“H(j«)J - 1 - (1 

1 + A6 


) f(co) 


(5. 60) 


where 


oo J < a 


<o| > O' 


(5. 61) 


A familiar way of approximating the rectangular function is with Butterworth 


functions of the form 


f(co) = 


n = 1.2, 


(5. 6 2) 


1 +( 5 > 


It turns out that we will achieve much better performance by replacing a by 
a parameter ft which we will vary to optimize the loop. With this change, and 
substitution of Eq. (5.62) into Eq. (5.60) we have 
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1 + A6 2 


(5. 63) 


!l - H(jco) | 2 


2n 

<!> + 


i + ( 


CO 

$ 


2n 
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This will approach the optimum given, by Eq. (5.58 ) both for very large and 
very small co, with the speed of convergence controlled by n. A particular 
advantage of our choice of Butterworth functions is that in calculating Ii(j<o) 

(and F(jco)) we can take advantage of the fact that the factors of terms of the 

2n ^ 

form 1 + x have been previously studied. 

Performance of these approximations has been calculated for n= 1, 2, 3 
and 4 and is shown in Fig. 16. These results are for a threshold defined by 
VAR [ e (t)] = 1/10. 

Another approaches to calculate the total loop error (again assuming 
rectangular input spectrum) for the second order loop discussed in section 5. 3. 
The loop parameters are then chosen to minimize the loop error variance. This 
gives the result shown with a dashed line in Fig. 16. The fact that these two 
“optimized" second order loops give essentially the same threshold suggests 
that we have roughly established the limit to performance obtainable with a 
second order loop. However, we can see from Fig. 16 that significant im- 
provements can be achieved with higher-order loops. For instance, with' a 
modulation index of 10 the fourth order loop gives 2 db of improvement over 
the second order loop and comes within l 1 ^ db of the optimum. 


See, for example, p. 252 of Seshu and Balabanian^ Ref. 3. 
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A. = Intrinsic Signal -to- Noise Ratio, (db) 


Fig. 16 Performance of approximations to the optimum transfer function. 

(Var [e £ (t)] = 0. 1) 
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VI. TECHNIQUES FOR EXPERIMENTAL MEASUREMENTS 

6- I Previous ADCOM Work on Signal -to -Noise Ratio Measurement 

Extensive laboratory study on the techniques of signal -to -noise ratio 

measurement has previously been performed at ADCOM. Results of this work 

1 2 

have been reported in ADCOM reports. * The earliest of these measurement 
systems, shown in Fig. 17 provided useful data on the operation of an oscil- 
lating limiter at 455 kc. 

Another system, shown in Fig. 18 was based on coherent separation 

of signal and noise components. The results obtained with this system were 

2 

not so consistent as the previous system and the approach was dropped. 

Further work on signal-t'o-noise ratio measurement is reported in Ref. 1. 
The functional diagram of the system employed is shown in Fig. 19. This 
system measures the r-f signal -to-noise ratio and is usable after nonlinear 
predetection signal-with-noise processing. The present measurements do not 
require the use of this complex system since all linear circuits may be employ- 
ed up to the detector circuits requiring measurement. 

The system designed for the present project and currently under con- 
struction relies on the techniques which have proven effective in these previous 
studies. In particular, the first design mentioned above provides much of the 
basis for the present design. 
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Fig. 17 Test system for the measurement of S/lST ratio. 
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Fig. 18 S/N meter and block diagram of the experimental configuration. 
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Fig. 19 Functional diagram for the ADCOM S/N meter 












The demodulator threshold test system is required to measure char- 
acteristics of the various FM demodulation schemes. The desire to obtain 
accurate results requires that the system allow for calibration repeatability 
to within one -half of a db {+ 5%) overall. To achieve this accuracy requires 
that the short-term (15-30 minute) stability of many of the individual com- 
ponents be + 0. 1 db or better. The long-term stability need not be ofhighac- 
curacy because provision is made to calibrate the system at intervals. However, 
the requirements for calibration should be held to a minimum by providing 
monitoring circuitry. 

One method of providing calibrated signal and noise levels is shown in 
the diagram of Fig. 20. The FM signal is generated, passed through a 
precision calibrated attenuator, and added to a noise source that has also been 
passed through an attenuator. The composite signal plus noise is filtered in 
stable measurable filters to provide precise noise bandwidths. 

The filtered composite signal is fed simultaneously to a wide band- 
width-limiter-discriminator (to be used as a standard of comparison for threshold 
improvements) and to the demodulator under test. The resultant demodulated 
signal is filtered with a calibrated lowpass filter. This provides a known base- 
band noise bandwidth for the precision measurements. The output signal level 
is measured with a narrowband spectrum analyzer and the noise level is meas- 
ured with a true rms meter after the signal has been removed with a notch 
filter. 
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{ Qucm-Tech 303 ) 



Fig. 20 Manual demodulator threshold test system 















Measurements on this system are taken by adjusting the input signal - 


to-noise ratio with the step attenuators in steps of one db. The output signal 
and noise levels are read and recorded. The threshold curve is plotted by 
calculating the signal -to -noise ratio at the output from the meter readings for 
each point on the curve. 

6. 2. 1 Automatic Threshold Curve Plotter 

The addition of certain laboratory equipment allows the time 
consuming point-by-point method to be replaced by a rapid and more accurate 
automatic .plotting system. This system is shown in Fig. 21. 

The basic system remains unchanged with two exceptions. The 
step attenuator in the noise channel is replaced by a motor driven attenuator 
and the signal portion of the demodulator output is held constant by an AGC 
circuit. Provision is made for applying a d-c voltage proportional to the in- 
put signal -to -noise ratio (in db) to the x input of an x-y plotter. The y input 
of the plotter is taken from the true rms noise measurement and passed through 
a d-c logarithmic amplifier so that it is proportional to the output signal -to -noise 
ratio in db. 

The required additions to the system in order to provide the 
greatly increased accuracy and repeatability of the automatic system are 

a) the motor driven attenuator, and 

b) the AGC control circuit. 

The other components of the system are standard commercial 
laboratory equipment. 
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Fig. 21 Automatic demodulator threshold test system. 
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VII. OUTLINE OF WORK FOR THE NEXT INTERVAL 


The work planned for the next Interval may be divided into a 
theoretical part and an experimental part. 

7. 1 Outline of Theoretical Work 


This part of the planned work, expected not to exceed- one third 
of the effort, will be as follows: 

a) Complete the evaluation of the work by Akima, Battail, 
Schilling and Billig, Van Trees and Viterbi. 

b) Complete the analysis of the band- dividing technique 
based on sample -and -hold message representation to 
determine 

(i) Effects of filter mismatch and detuning, and 
the associated AM and FM transients giving 
rise to inter -sample interference and erron- 
eous estimation of the frequency. 

(ii) Devise and refine schemes for minimizing 
effects in (i). 

(iii) Determine the effects of phase steps in going 
from one filter output to the other (the in- 
coherency problem) and evaluate methods for 
avoiding these steps or abating their effects. 

c) Complete the analysis of combined PLL and FCF de- 
modulation technique. 

d) Complete evaluation of other approaches specified in 
the Work Statement. 

7. 2 Outline of Experimental Work 

This part of the planned work, expected to amount to a major 
fraction of the effort, will be as follows: 

a) Complete laboratory testing preparations • 
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b) Complete design evaluation analyses of the new 
implementation concepts, developed on this 
program, of the band- dividing technique and 

• related circuitry. 

c) Complete design and implementation of a 
higher -order phase -lock demodulator for early 
preliminary tests, and for providing a test bed 
for other promising variations on the basic 
phase -lock technique. 
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